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Abstract 

We consider a massless scalar field in 1+1 dimensions that satisfies a Robin 
boundary condition at a non-relativistic moving boundary. Using the perturbative 
approach introduced by Ford and Vilenkin, we compute the total force on the mov- 
ing boundary. In contrast to what happens for the Dirichlet and Neumann bound- 
ary conditions, in addition to a dissipative part, the force acquires also a dispersive 
one. Further, we also show that with an appropriate choice for the mechanical fre- 
quency of the moving boundary it is possible to turn off the vacuum dissipation 
almost completely. 

The interaction between a physical system and a material plate (or cavity in general) 
in its surroundings has a long history. In 1948, Casimir and Polder 1 1 1 computed for 
the first time the retarded interaction energy between a neutral but polarizable atom 
and a perfectly conducting wall. At this same year, Casimir 1 2 1 predicted the attraction 
between two neutral parallel conducting plates due to the shift caused by the plates in 
the energy of the radiation field in vacuum state. Casimir's result may be considered 
the first problem worked out in detail of the so called cavity QED. Since then, a lot of 
work has been done on the Casimir effect, see for instance the reviews I3ll4ll5ll6ll71l8l 
and references therein (for other phenomena of cavity QED, see I9l ll0h . 

However, the interaction between a quantum field and a material plate is quite com- 
plicated. Hence, as a first approximation, it is common to simulate this interaction by 
imposing an idealized boundary condition on the field . The most familiar conditions 
are Dirichlet and Neumann ones. A less familiar, but not less important condition is 
the so called Robin boundary condition, defined for a scalar field by 

deb 

<j>\dn = Pg-\dK, (1) 

where dlZ is the boundary of the system under study, |^ means n • V(f>, with n being 
a unitary vector normal to the boundary and j3 is a parameter with dimension of length 
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that can assume any value in the interval [0, oo). Robin BC have the nice property of 
interpolating continuously Dirichlet and Neumann ones. From 0, we immediatly see 
that for j3 — we have Dirichlet BC and for (3 — > oo we have Neumann BC. 

In this work, we discuss some consequences of using Robin BC in the context of 
the Dynamical Casimir effect. However, before starting our calculations, we shall make 
a few comments about this kind of BC. Robin BC already appear in a natural way in 
classical physics. For instance, when we solve problems in classical electromagnetism 
in the presence of spherical conducting shells the radial functions satisfy Robin BC 
with particular values of parameter j3. Another nice example, still in the context of 
classical physics, is the problem of a vibrating string subjected to a tension T with two 
massless rings at its ends which may slide without friction along vertical rods and are 
coupled to springs of constants Kiand k%, respectively, as indicated in figure^ 




Figure 1 : Elastic supports at x = and x = a give rise to Robin BC. 



Assuming small inclinations << l), application of Newton's Second Law to 
both massless rings gives 



T dy 
y\x=o — — tt- x=o 

Kl OX 



> I T dy 
and y\ x =a = -^-\x=a 

K2 OX 



(2) 



The fact that Robin BC simulates an elastic support at the boundary has been 
pointed out in the literature 11 II . Though the reflection at a fixed boundary where 
the wave satisfies a Robin BC is complete, there is some kind of time delay caused by 
a bulk/boundary dynamics. In other words, the reflection coefficient can be written as 



7? 



J4(k) 



(note that \R\ = 1), where k is the wavenumber of the incident wave and 



hence there will be a phase shift between the incident and reflected waves. This gives 
a qualitative explanation for the surface terms that appear in connection with Robin 
BC in quantum field theory H2lll3lfT4lll5l . Total energy (string plus surface terms) 
is conserved, but there is a "bulk/boundary" exchange, so that the energy of the string 
itself is not conserved: 



d_ 

dt Jo 



1 fdy\ 2 1 fdy^ 



dx 



«iy(0,t)^(0,t) 



, . dy , . 
K2y{a,t) — {a,t) 



Robin BC are also useful for phenomenological models that describe penetrable 
surfaces 1191 . In fact, for some particular cases, these conditions can simulate the 
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plasma model for real metals. It is not difficult to show that for frequencies much 
smaller than the plasma frequency, uj << up, a small value of ft plays the role of 
1 / ujp (c = 1). In other words, under such assumptions, ft is proportional to the plasma 
wavelength, which is directly related to the penetration depth of the field. 

Recently, Robin BC have been studied in many different contexts, namely: Bon- 
durant and Fulling 1161 discussed in detail the Green's functions of the wave, heat and 
Schrodinger equations under Robin BC; Albuquerque and Cavalcanti 1121 and Albu- 
querque |:243 analized the one-loop renormalization of a A(/> 4 theory under these con- 
ditions; Minces and Rivelles used them in the context of AdS/CFT correspondence; 
Solodukihn 1201 studied upper bounds for the ratio between entropy and energy of sys- 
tems constrained by Robin BC; heat kernel coefficients were studied by Bordag et al 
1221 . Fulling 1131 and Dowker 1231 : and a very detailed calculation of the static Casimir 
effect with Robin BC was made by Romeo and Saharian 1151 . It is worth mentioning 
that Robin BC may give rise to restauring Casimir forces between two parallel plates, 
once parameters ft at each plate are appropriately chosen. 

However, since the pioneering paper by Moore |25'] on radiation reaction forces on 
moving boundaries, Robin BC have never been considered explicitly in the context of 
the dynamical Casimir effect (as far as the authors know). It is our purpose here to 
make this kind of calculation in a simple model, namely, we shall consider a massless 
scalar field in 1+1 dimensions subjected to a Robin BC at one non-relativistic moving 
boundary. The main motivation is the following: it has been shown that for Dirichlet 
1261 1271 1281 and Neumann BC 1291 the linear susceptibilities are equal and purely 
imaginary, 

x ° { u>)= X »{u)=i^. ( 3 ) 
D7rcr 

These susceptibilities lead to purely dissipative forces on the moving boundary: 

5F°( t ) = 5F»(t) = J^P q (t), (4) 

where 5q(t) is the position of the moving boundary at instant t. 

For more general BC see Jaekel and Reynauld 1281 1301 and for 3+1 calculations 
see references 1311 . Since Robin BC interpolates continuously Dirichlet and Neumann 
ones we are led to make the following questions: what happens to the force for the 
interpolating BC? Will it still be a purely dissipative one? In what follows we shall 
answer these questions. 

Besides the assumption of a non-relativistic motion for the boundary, we shall also 
suppose that the boundary has a prescribed motion with a small amplitude, 5q(t) being 
its position at time t. Hence, we assume that 

\6q(t)\ < c and \Sq(t)\ < c/ujq , (5) 

where ujq corresponds to the typical mechanical frequency. Therefore, we need to solve 
the wave equation for que quantum field, d 2 (f>(t, x) = 0, with <fr satisfying a Robin BC 
at the moving boundary, which, in the co-moving frame, is written as 

-^j(t' ,x')\Bound. = -p(j>'(t',x')\Bound. ■ (6) 
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The corresponding BC in the laboratory frame is given by 



9 "3 
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<t>{t,x)\ x =s g ( t ) = -<j>(t,x)\ x =s g ( t ) , 



(7) 



where we neglected terms of 0(8q 2 / c 2 )). Using the perturbative approach of Ford and 
Vilenkin 1 27 1 we write 



cf>{x,t) = <j> {x, t) + 5(t>{x,t) , (8) 
where </>o is the solution with a static boundary at x = 0, which is given by 

<Po{t,x) 



sm(u) x) + uif3 cos(uix) a(ui)e %ut + a) (u>)e 



(9) 

and 5cf) corresponds to the contribution generated by the movement of the boundary. 
This perturbation satisfies the wave equation d 2 5(f>(x, t) = with the following BC: 



-&K*)-^r(*,0)> (10) 



where we discarded terms of 0(<5q 2 ). The total force on the boundary is given by: 

5F(t) = (0\T n (t,Sq + (t)) — T u (t,5q~(t))\0) , (11) 

where T n (t,x) = -|i (d x (f)) 2 (t,x) + (d t (t>) 2 (t,x) I . Substituting = 4> + 5<jr. 



SF(t) = -^(0\(^{(d x cf >0 )(tJq + (t)),(d x 5cf > )(t 1 Sq+(t))} + 

+ { (dtfo) (t, 8q+ (*)), (9 t ty) (t, $g+ (<))}- [5q+ (t)^5q-{t)]\ |0) + 

In the last equation {•■,■■} means anticomutator and terms involving only the non- 
perturbed field (f>o disappear. Now, we expand around x = and keep only first order 
terms. One may also show that the total force is twice the force on each side. With 
these facts in mind, we get 

5F(t) = -i(O|^{(9^o)(i,O + ),(a^0)(t,O+)} + 

+ {(a t </> )(t,o + ),(a t ^)(t,o+)}- [o+-o-]j|o) 
= -(o|{(0 x 0o)(t, o + ), (dM){t, o+)} + {(at0o)(t, o + ), (W)(t, o+)} |o) , 

Denoting by 8F(lo), 6&(cj,x) and5Q(uj) the time Fourier transforms of SF(t), S(j)(t,x) 
and <5g(i), respectively, it is straightforward to show that 



6F(lo) 



~2 



-I (0\[d x $ o {LJ-u;',0),d x S$(w',0)} 
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- (w-L)')w'{<?>o(w-v',0),S$(u',0)}j |0) . 

Hence, we must solve the equation (d 2 + oj 2 )8$(x, u) = with the BC (this is condi- 
tion translated to the Fourier space): 

1 1 r J, . // 

'u "$ (^ ", 0)SQ(u '-u"). 

Z7T 

However, <5<£>(w, x) satisfies a second order differential equation, which means that we 
shall need an extra condition. A natural choice is to consider only the solutions for 
64>(t, x) which describe perturbations getting away from the boundary: 

6$(u,x) = 8gn(x)—d x 6$(uj,0)e iu W => 5$(uj, ± )= ±— d x 6$(cj, 0) . 

(12) 

Last equations allow us to express 8$(u>, 0) and d x 8$>(u>, 0) in terms of the static field. 
The resulting expressions, when substituted in 8J-{uS), give 

fdu'f J3u/ \ fduj" , 
- i(0|{5 x $o(w-w',0),a x $o(w",0)}|0)- 

+ ... + (w-w / )(w / -w"V"(0|{$o(w-cj' ! 0),$o(w",0)}|0)V 

where ... means other (though analogous) correlators. However, all correlators in the 
previous expression are connected by the field equation and Robin BC to the follow- 
ing one: Cbfwi^a) = (0|{$o(wi, 0), $o( w 2, 0)}|0), which involves only the non- 
perturbed field. A straightforward calculation leads to 



C {uji,uj 2 )= ( -] , 2 2 . |cji|(5(a;i +cl> 2 ) ■ (13) 



4tt/3 2 
(1 T^f/3 2 ) 

With the aid of this correlator, we write 8T(w) in the form 

&FM =: <*QM , d4) 
where the real and imaginary parts of the susceptibility x( w ) are identified as 

^ xM = [^'-^+1]^+!) (15) 

and 

1 oj'\uj'-uj\(i + 2[3 2 u;'(uj' -uj) + (3 4 l;' 2 (u J ' -uj) 2 } 

JtoyM=-/ du/ i = n -. (16) 

1 ' kJ-oo [/3 2 (^'- W ) 2 + 1](/3 2 o;' 2 + 1) 
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Before we proceed, it is interesting to check some limits. Taking the limits /3 = 
(Dirichlet BC) or j3 — ► oo (Neumann BC) in the above expressions, we obtain 

i f +oc 

lZex{u) — > and lmx{u) — ► — duj ' oj ' \oj — oj ' \ . (17) 

T ./-oo 

As anticipated, the susceptibility is purely imaginary in these limits. In order to perform 
the integration for Xmx(w) we need a regularization prescrition. In this case, a very 
natural way of doing that is to write the integral in the form 



dw'w'\w-w'\\ = Jim [ / + + \<L)'u'\u-u'\. (18) 

—L JO Jlu ) 



L — >oo 



Note that the first and third integral on the right hand side of the above equation cancel 
out. Therefore, we are left with 

h c u hoi 3 
X(w) = i-/ du'u'\u-u , \=i— , (19) 

which leads to the well known forces already written in (0}. For later convenience, it is 
worth emphasizing that the total work is given by 

/+oo rCO 
F(t)Sq(t)dt= / dLOLulm X (uj)\SQ(uj)\ 2 . (20) 
-oo 7T Jo 

Note that, for Dirichlet and Neumann BC, Imx(w) > for to > 0, so that the force is 
always dissipative. 

Using a regularization prescription analogous to that described above in equations 
dl5> and d!6i . the contributions for the integrals coming from the intervals (— oo, 0) and 
(ui, oo) will cancel each other and we are left with integrals from to ui. Performing 
the remaining integrals, we obtain: 

„ , s u -(3Lu\og<P 2 cu 2 + I) + I3 3 cu 3 + A/3lu -2lg- 1 ( (3lu){P 2 lu 2 + 2) 

Ue X(u) = — 2 R2 2 i A 

p z 7T p z LU z + 4 

^ . , u /3 4 w 4 + 4/3 2 w 2 -61og(/3 2 w 2 + l)( / 3 2 w 2 + 2) + 12/3cjtg- 1 (/?u;) 
lm X (w) = —2 tt2 2 T a ■ 

Expanding the previous expressions appropriately, the first corrections to the Dirichlet 
and Neumann cases can be obtained: 

For /3u> << 1, we have 

£.4 2ll^ uj 3 uj^ 

He X (uj) = -—[3+—P 3 + 0({3 5 ) and Im X {u) = — - — (3 2 + 0((3 4 ) 

07T l07T D7T D7T 

For (3lu >> 1, we have 
ne X (.) = -^l-^ + 0(^) and lm X (.) = £ - H^^M + O^) 

7T [j j.V D7T TT [V 
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Hence, the total force is given by 

SF(t) = ±^5q(t)-f3^5 q (t)}+0(f3 2 ) (/J -0) 

The behaviours of 1Zex(u) and Imx(w) are shown in figure D 
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Figure 2: Real and imaginary parts of x{lo), conveniently normalized by Imxc(w), 
as functions of f3u. 

In this work, we computed the total force on a non-relativistic moving bound- 
ary in 1+1 dimensions due to the vacuum fluctuations of a massless scalar field sub- 
jected to Robin BC. Dirichlet and Neumann BC correspond to particular limits of our 
results (particular values of /?). It is worth emphasizing that when Robin BC are 
used the susceptibility acquires a real part. The pronounced valley in the graph of 
Imx(w)/Tmx£)(w) leads to a quite interesting result: if SQ(u>) is peaked around ujq, 
equation J20i shows that, for any fixed (3, there will be an appropriate choice of ojq 
such that the dissipative effects on the boundary will be almost completely eliminated. 
A natural sequence of this work is to compute the particle creation rate under the same 
circumstances as those assumed in this work. This problem is under study and the 
results will be published elsewhere. 

The authors thank CNPq and FAPERJ for financial support and S.A. Fulling for 
bringing reference 11 II to our knowledge. 
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